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ABSTRACT 

We present the discovery of a transiting exoplanet candidate in the K2 Field-1 with an orbital period of 
9.1457 hr; K2-22b. The highly variable transit depths, ranging from ~0% to 1.3%, are suggestive of a planet 
that is disintegrating via the emission of dusty effluents. We characterize the host star as an M-dwarf with 
Teff ~ 3800 K. We have obtained ground-based transit measurements with several 1-m class telescopes and 
with the GTC. These observations (1) improve the transit ephemeris; (2) confirm the variable nature of the 
transit depths; (3) indicate variations in the transit shapes; and (4) demonstrate clearly that at least on one 
occasion the transit depths were significantly wavelength dependent. The latter three effects tend to indicate 
extinction of starlight by dust rather than by any combination of solid bodies. The K2 observations yield a 
folded light curve with lower time resolution but with substantially better statistical precision compared with 
the ground-based observations. We detect a significant “bump” just after the transit egress, and a less significant 
bump Just prior to transit ingress. We interpret these bumps in the context of a planet that is not only likely 
streaming a dust tail behind it, but also has a more prominent leading dust trail that precedes it. This effect 
is modeled in terms of dust grains that can escape to beyond the planet’s Hill sphere and effectively undergo 
‘Roche lobe overflow,’ even though the planet’s surface is likely underfilling its Roche lobe by a factor of 2. 
Subject headings: planetary systems—planets and satellites: detection, atmospheres 
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1. INTRODUCTION 

The Kepler mission (Bomcki et al. 2010) has revolutionized 
the field of exoplanets, with some 4000 planet candidates dis¬ 
covered to date (Mullally et al. 2015), of which at least 1000 
have been confirmed (Lissauer et al. 2014, Rowe et al. 2014). 
With the original objective of discovering Earth-size planets 
in the habitable zone of their host stars, the telescope was 
bound to also improve our understanding of close-in rocky 
planets (Jackson et al. 2009; Schlaufman et al. 2010). In¬ 
deed, the first Kepler rocky planet, Kepler-10b, had an orbital 
period of only 20 hr (Batalha et al. 2011). The smallest planet 
with a well measured mass and radius, Kepler-78b, also has a 
very short orbital period of 8.5 hr (Sanchis-Ojeda et al. 2013), 
which was instrumental in measuring its small mass of 1.7 
Earth masses (Howard et al. 2013, Pepe et al. 2013). In spite 
of the falloff in the numbers of Kepler exoplanet candidates 
at short periods, there are 106 well vetted candidates with or¬ 
bital periods shorter than one day (hereafter “USPs”; Sanchis- 
Ojeda et al. 2014), and most of them seem to be smaller than 
twice the size of Earth. 

Not included among the above lists are two special transit¬ 
ing exoplanets that are thought to be disintegrating via dusty 
effluents (Rappaport et al. 2012; Rappaport et al. 2014). In 
both cases it is inferred that the planets are trailed by a dust 
tail whose dynamics are influenced by radiation pressure on 
the dust grains. This leads to transit profiles characterized by 
a pronounced depression in flux after the planet has moved 
off of the stellar disk (i.e., a post-transit depression). In the 
case of KIC 12557548b (Rapapport et al. 2012; hereafter 
‘KIC 1255b’) the transit depths range from ~1.2% down to 
<0.1% in an highly erratic manner, while for KOI 2700b 
(KIC 8639908; Rappaport et al. 2014) the transit depths are 
observed to be slowly decreasing in depth over the course of 
the fours years of Kepler observations. The fact that these 
‘disintegrating’ planets are relatively rare (2 of 4000 Kepler 
planets) is likely due to the conditions required for their ex¬ 
istence and detection, namely high surface equilibrium tem¬ 
peratures and very low surface gravity, and a possibly short 
disintegration lifetime of only 10-100 Myr (see, e.g., Rappa¬ 
port et al. 2012; Perez-Becker & Chiang 2013). 

The main Kepler mission had an abrupt ending when two 
reactions wheels failed by March 2013. The reaction wheels 
are very important to maintain the telescope pointing in a 
given direction, and the telescope could no longer point to¬ 
ward the original Kepler field. The problem was partially by¬ 
passed by designing a new mission, called “K2”, in which the 
telescope would point toward a different field of view along 
the ecliptic plane every three months (Howell et al. 2014); 
the spacecraft stability is improved by equalizing the Sun’s 
radiation pressure forces on the solar panels. The unfortunate 
demise of the reaction wheels that put an end to the main mis¬ 
sion, also opened the possibility for new discoveries of planets 
orbiting brighter stars since thousands of new bright stars are 
observed in each field. 

In the short lifespan of this new mission, there have been 
several papers describing techniques to produce light curves 
(Vanderburg & Johnson 2014; Aigrain et al. 2015; Eoreman- 
Mackey et al. 2015), and planet discoveries like a super- 
Earth transiting a bright host star (HIP 116454, Vanderburg 
et al. 2015), a triple planet system orbiting a bright M-dwarf 
(K2-3, Crossfield et al. 2015), and a pair of gas giants near a 
3;2 mean motion resonance (EPIC 201505350, Armstrong et 
al. 2015b), with almost 20 confirmed K2 planets discovered 


to date (Montet et al. 2015). There have also been several 
catalogs of variable stars and eclipsing binaries (Armstrong 
et al. 2015a; LaCourse et al. 2015). This paper is the first in 
a series describing our discoveries using the K2 public data 
releases. The name of the project, “ESPRINT”, stands for 
"Equipo de Seguimiento de Planetas Rocosos INterpretando 
sus Transitos", which in English means "Eollow-up team of 
rocky planets via the interpretation of their transits". 

In this work we focus on the surprising discovery of another 
one of these candidate disintegrating planets, this one in the 
K2 Eield 1 which contains only 21,647 target stars (close to 
an order of magnitude fewer than in the prime Kepler field). 
Even more impressive, this particular short-period exoplanet 
appears to have a dominant leading dust tail and possibly an 
additional trailing one, a phenomenon not seen before in as¬ 
trophysics. The paper is organized as follows. In Section |2] 
we summarize the observations taken with the K2 mission 
and describe how this particular object was found. In Sec- 
tion[3]we describe the variable transit depths, the timing anal¬ 
ysis, and the unusual transit profile that cannot be explained 
by a solid body. We present and discuss 15 transit measure¬ 
ments that were made in follow-up ground-based observations 
in Section m We analyze the properties of the host star based 
on a number of ground-based imaging and spectral observa¬ 
tions in Section |5] In Section |6] we set significant constraints 
on the radial velocity variations in the host star. We discuss 
the wavelength dependence of the transits observed with the 
GTC in Section|2] We summarize why the host of the transits 
is the bright target star and not its much fainter companion 
in Section [8] In Section |9] we interpret all the observations 
in terms of a model in which the planet is disintegrating, and 
discuss why the different characteristics and environment of 
K2-22b could lead to a dominant leading dust tail. Einally, in 
Section [To] we present a summary and conclusions and point 
toward new lines of research that could improve our under¬ 
standing of how disintegrating planets form and evolve. 

2. K2 DATA PROCESSING 

The target star, with EPIC number 201637175 (from now 
on named K2-22), was selected as one of the 21,647 stars in 
Eield 1 to be observed in the long cadence mode of the K2 
mission (Howell et al. 2014). During the period from 2014 
May 30 to 2014 August 20, a total of 3877 images of 15 x 15 
pixels were recorded by the Kepler telescope, with a typi¬ 
cal cadence of 29.42 minutes. The data were sent to NASA 
Ames, subsequently calibrated, including cosmic ray removal 
(Howell et al. 2014), and uploaded to the public K2 MAST 
archive in late 2014 December. The data were then down¬ 
loaded from the MAST archive and utilized for the analysis 
presented in this work. 

The discovery of K2-22b is part of the larger ESPRINT col¬ 
laboration to detect and quickly characterize interesting plan¬ 
etary systems discovered using the K2 public data. In this sec¬ 
tion we highlight the way in which we produce light curves for 
all the observed stars, how this object was identified as part of 
the survey, and how we produced a better quality light curve 
for this particular object once the transits had been detected. 

2.1. ESPRINT Photometric Pipeline 

Our photometric pipeline follows the steps of similar ef¬ 
forts published to date (Vanderburg & Johnson 2014; Cross¬ 
field et al. 2015; Aigrain et al. 2015; Eoreman-Mackey et al. 
2015; Lund et al. 2015) that describe how to efficiently ex¬ 
tract light curves from the calibrated pixel level data archived 
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on MAST. The ingredients to generate the light curves, with 
our own choice described, are: 

• Aperture selection: Our apertures have irregular 
shapes, which are based on the amount of light that 
a certain pixel receives above the background level. 
We selected this type of aperture to capture as much 
light as possible while reducing the number of pixels 
used, which in turn reduces the noise induced by a large 
background correction. Based on experiments we car¬ 
ried out on the engineering data release, selected pixels 
must be 30% higher than the background level (esti¬ 
mated from an outlier corrected median of all the pix¬ 
els in the image) in 50% of the images in the case of 
a star brighter than Kepler magnitude Kp = 11.5. For 
stars fainter than magnitude Kp = \A, the selected pix¬ 
els must be 4% higher than the background. For stars of 
intermediate magnitudes, a linear interpolation of these 
two thresholds is used. A simple algorithm groups con¬ 
tiguous pixels in different apertures , and the target star 
aperture is selected to be the one that contains the tar¬ 
get star pixel position (obtained from the FITS head¬ 
ers). This type of aperture is similar to the ones used by 
Lund et al. (2015), and quite different from the circular 
apertures used in many of the other pipelines. 

• Thruster event removal: As highlighted in Vanderburg 

6 Johnson (2014), every 6 hours the telescope rolls to 
maintain the targets on the defined set of pixels that are 
downloaded, in what is known as a ‘thruster event’. 
We recognized these events by calculating the cen¬ 
troid motion of a particularly well behaved star (EPIC 
201918073), and selected those moments where the po¬ 
sition of the star jumps much more than usual (in the 
case of these stars, 0.1 pixels in the x direction). The 
images obtained during thruster events are removed 
from the analysis. 

• Data slicing: We split our dataset into eleven differ¬ 
ent segments chosen to have a length of approximately 

7 days, but also to contain an integer number of tele¬ 
scope roll cycles. The first segment and the one after 
a large gap (in the middle of the dataset) are not used 
in the global search since they are poorly behaved in 
some cases, with systematic effects induced by thermal 
changes that appear after reorienting the telescope. 

• Systematics removal: We calculate the centroid po¬ 
sitions and obtain a fourth-order polynomial that de¬ 
scribes the movement of the star in the x and y coordi¬ 
nates. This polynomial fit is used to determine a new set 
of coordinates, in which the star moves only along one 
direction. The fluxes are then decorrelated first against 
time and then against this moving coordinate, using a 
fourth-order polynomial in each case. This process is 
repeated 3 times, and the results are very robust against 
problems caused by the presence of low-frequency as- 
trophysical sources of noise (see Vanderburg &. John¬ 
son 2014 for a more detailed description on how this 
process works). 

This recipe was followed to generate the light curves of the 
21,647 stars in Field-1. Among them, the light curve for EPIC 
20163717 can be seen in the top panel of Pig. [T] It is interest¬ 
ing to note that since our apertures depend on the amount of 


background light, which is increasing over the course of the 
observations, the number of pixels in the aperture decreases 
with time. This is a desired effect, since it tends to balance 
the increase of background light in the aperture by reducing 
the number of pixels, and therefore changes in the flux scatter 
are less severe. 

These light curves are generally analyzed using two differ¬ 
ent search algorithms: a more standard BLS routine (Kovacs 
et al. 2002; Jenkins et al. 2010; Oflr 2014) to search for plan¬ 
ets with orbital periods longer than 1 day, and a more special¬ 
ized PET pipeline used to detect planets with orbital periods 
shorter than l-2days (Sanchis-Ojedaetal. 2014). In this case, 
due to the short orbital period of the signal, we describe only 
the PET search. 

2.2. Detection via the FFT technique 

Our target star K2-22 was detected as part of our search 
for ultra-short period planets in the K2 Pield-1 dataset using 
the FPT technique. The routine automatically identifies those 
objects for which a main frequency and at least one harmonic 
can be distinguished above the level noise in the PET power 
spectrum (see Sanchis-Ojeda et al. 2014 for details). A total 
of 2628 objects were identified in that way, but a large fraction 
of them were caused by improper corrections of the 6-hour 
roll of the Kepler telescope. These false detections are easy 
to remove since their main frequency is always related to the 
fundamental roll frequency of 4.08 rolls per day, although this 
simple removal clearly affects the completeness of our search. 
A total of 390 objects were selected for visual inspection, and 
among them K2-22 was selected as the most promising ultra- 
short orbital period planet candidate. 

The PET by which this object was discovered is shown in 
Fig.E] Note the prominent peak at 2.62 cycles/day which is 
the base frequency corresponding to the 9.1457-hour period, 
as well as the next 8 higher harmonics which lie below the 
Nyquist limit. The overall slowly decaying Pourier ampli¬ 
tudes with harmonic number is characteristic of short-period 
planet transits. 

2.3. Individualized aperture photometry 

Our method for generating the light curves relies on a one- 
to-one relationship between the raw flux counts and the posi¬ 
tion of the star on the CCD chip. Any source of astrophysical 
variability could distort this relationship, and this is the case 
for both stellar activity induced signals and transits. A closer 
inspection of the raw light curve of K2-22 shows long-term 
trends that do not correlate with the centroid motion, and are 
likely due to the slow rotation of the host star. These trends 
are removed automatically as we fit a fourth-order polynomial 
in time to each of the 7-day segments, which effectively re¬ 
moves variability on scales longer than approximately 2 days 
(see upper panel of Pig. [Til. 

In order to remove the effects of the transits, we first folded 
the original light curve given the period obtained from the 
PFT, after which we identified those orbital phases where no 
transit is expected. We then ran our photometric pipeline 
again, using the same aperture (see Pig.|3]l but only using the 
out-of-transit flux measurements to And the best fit polynomi¬ 
als to correct for both temporal and telescope motion varia¬ 
tions. This process reduced the photometric scatter, and en¬ 
couraged us to try different approaches to continue improving 
the light curve. We tried different combinations of polynomial 
orders and also different approaches to defining the apertures. 
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Fig. 1. — Lightcurve of K2-22 in time bins equal to the half-hour Kepler 
long-cadence sampling time. Top panel: Light curve of the form used in our 
global search for USP planets. Middle panel: Light curve processed with a 
modified algorithm that better preserves stellar activity. Lower panel: Au¬ 
tocorrelation function, with a vertical line representing the infeiTed rotation 
period of the star. 


but none of them improved the quality of the light curve (see 
top panel of Fig.[T]i. After all the corrections, the final scatter 
per 30 minute cadence is 650 ppm, which is near the mean un¬ 
certainty obtained with our photometric pipeline for the typi¬ 
cal K2 15* magnitude star. 

We also tried to produce a light curve in which other astro- 
physical signals would be preserved (e.g., starspot rotation). 
During the process of detrending each of the 11 segments of 
data, we saved the coefficients of the fourth-order polynomial 
in time, and used them to reconstruct the signal again after 
removing the centroid motion artifact. Since the aperture is 



5 10 15 20 

Frequency (cycles/day) 

Fig. 2.— Discovery Fourier transform of the flux data showing strong peaks 
at the 9.147-hour period and all 8 hai'monics that are below the Nyquist limit¬ 
ing frequency. The arrows mark two additional harmonics aliased around the 
Nyquist limit. 



Fig. 3.— K2 image of the object K2-22. A customized aperture is defined 
based on the amount of light of each pixel, and level of background light. The 
blue star represents the expected position of the target star, and the Kepler 
magnitude obtained from the EPIC catalog. 

individually defined for each segment, we had to adjust the 
mean flux level of each segment to create a continuous light 
curve. This astrophysically more accurate light curve is also 
shown in Fig.[T] and exhibits a clear signal of starspots with a 
rotation that could either be 7-8 days or twice this value. The 
shorter quasi-periodicity would typically arise when the star 
has two active longitudes separated by 180 degrees in lon¬ 
gitude. We used an autocorrelation function to confirm this 
suspicion (see lower panel of Fig. [U, and measured a rota¬ 
tion period of 15.3 days, following the techniques described 
in McQuillan et al. (2013). 

3. TRANSIT PARAMETERS AND DEPTH VARIABILITY 

In this section we describe the transits of K2-22b as ob¬ 
served by K2, with particular emphasis on the characteristics 
that deviate from the transits of a more typical planet. 
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Fig. 4.— Zoomed version of Fig. ^where the variability of the transit 
depths in much easier to see. Each vertical red line represents an expected 
transit time. 


3.1. Individual transit times and depths 

The top panel of Fig. [T] shows the full data set from the 
K2 Field-1 observations covering an interval of ^80 days. It 
is apparent that there are sharp dips in intensity whose depths 
are highly and erratically variable. In Fig.|4]we can see a zoom 
in on two weeks of observations. The individual transits are 
now quite apparent, and the depth variations are dramatically 
evident. 

In order to analyze these variations, we first folded the light 
curve with the period obtained from the FFT, after removing 
any long-period signals. We did this by fitting for local lin¬ 
ear trends using a total of 5 hours of observations right before 
and after each transit. We then fit the folded light curve with 
a simple idealized transit profile comprised of 3 straight-line 
segments: a flat bottom and sloping ingress and egress with 
the same slope magnitudes, hereafter referred to as the “three- 
segment (symmetric) model”. This is very similar to the sim¬ 
pler ‘box model’, but with non-zero ingress and egress times. 
We forced the fit to have a duration of at least 5 minutes on 
the lowest part of the transit, to make sure that the final fit did 
not look like a triangle. This fit gave a mean transit depth of 
0.5%, and a total duration of 1 hour and 10 minutes. Given 
that the cadence of the observations is approximately 30 min¬ 
utes, the real mean duration of the K2 transits must be close 
to 40 minutes, in agreement with what is expected of an USP 
planet orbiting an M-dwarf. This is confirmed by our follow¬ 
up ground-based observations with their better temporal reso¬ 
lution (see Sect. El. 

With the mean transit profile and a good estimate for the or¬ 
bital ephemeris, we fit each of the 190 individual K2 transits 
with the same three-segment model, allowing only the tran¬ 
sit time to vary. In the process, we evaluated the robustness 
of the detection of each of the individual transits. In some 
cases, the depth is so low that the transit cannot be detected, 
whereas on other occasions the transit occurs during a thruster 
event, so there are no data. In rare cases, the transit consists 
of only a single flux point, which does not allow for a clear 
determination of the transit time or depth. After removing all 


those cases, we were left with a reduced sample of 60 well 
detected transits. In addition to the best fit transit time, we es¬ 
timated the formal uncertainty by finding the interval of times 
where the standard function is within 1 of its minimum 
value, where a constant value of 650 ppm is used for the un¬ 
certainties in the flux. The transit times (including uncertain¬ 
ties) obtained in this manner are fit to a linear function, and 
the 0-C (‘observed minus calculated’) residuals of that fit are 
displayed in Fig.|5] A clear excess of scatter (above the formal 
statistical uncertainties) is detected, with a best-fit standard x^ 
of 480 for 60 transit times. Formal uncertainties, however, un¬ 
derestimate the true uncertainties when the model parameters 
are correlated, as could be the case here since we have re¬ 
moved a local linear trend for each transit which is known to 
be correlated with the transit time. We further examine this 
excess scatter below with simulations of the data train. 

We now use the new orbital ephemeris to fix the time of 
each transit, and repeat the process but now fitting for the 
transit depths. These transit depths are significantly differ¬ 
ent from one transit to the next, and the depths range from a 
maximum of 1.3% to 0.27% which is close to the photometric 
detectability limit of our time series for individual flux points. 
We can also see in Fig.|5]that the transit depths do not strongly 
correlate with the 0 — C timing residuals. These erratic tran¬ 
sit depth variations are quite reminiscent of those exhibited 
by KIC 1255b (Rappaport et al. 2012; Croll et al. 2014), and 
cover much the same range in depths. 

We checked for periodicities in the measured transit depths 
and timing residuals using a Lomb-Scargle periodogram, but 
no significant peak was detected with a false alarm probability 
smaller than 1%. We also attempted to constrain the true size 
of the planet by studying the shallowest K2 transits. We revis¬ 
ited the discarded transits, and selected the 6 shallowest cases 
in which the transit observations were complete (no thruster 
events). The mean depth of these transits is 0.14 ifc 0.03%, 
which, given the radius of the star (see section fOl i. translates 
into an upper bound of 2.5 ± 0.4 on the planet radius. 

Even though the transit depth variations described above 
are later confirmed by ground-based observations, there is a 
potential concern that the variability could have been caused 
by the relatively comparable transit and sampling timescales. 
In this same regard, it is also possible that some or all of the 
excess variations in the 0-C scatter (see Fig.|5]l over those 
expected from statistical fluctuations are due to the relatively 
short transit duration compared with the LC integration time. 
We have therefore carried out extensive numerical simulations 
of these effects. 

The simulations of the depth and transit-time variations are 
based on a model that has a simple box transit profile which 
is integrated over the 30 minute LC time, and includes a sinu¬ 
soid (and its first harmonic) to represent the rotating starspot 
activity. The amplitude of the sinusoids and the rotation pe¬ 
riod are fixed at 1% and 15 days, respectively. The duration of 
the model transit is fixed at 50 min (see Sect. E. The model 
is evaluated at the same times as the K2 observations, and 
only the same 60 transit windows are analyzed to ensure that 
the simulations represent a similar dataset to the one used in 
this paper. There are a total of four different numerical ex¬ 
periments, each one repeated many times, in which we either 
include starspots or not, and we either have a constant transit 
depth or depths drawn from a Gaussian distribution of a given 
variance, but always with a mean of 0.5%. In all cases a fixed 
orbital ephemeris is used. White noise of 700 ppm per LC 
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Fig. 5.— In depth analysis of individual K2 transits. Upper left panel: Residuals of the K2 transit times with respect to a linear ephemeris, where the error bars 
reflect only the formal uncertainties. The best solution has a final of 480 for 60 transit times. Our simulations show that this excess of scatter comes from 
systematic uncertainties induced by the poor time resolution compared to the duration of the K2 transits. Lower left: The individual transit depths with time; 
these confirm the erratic variations in the transit depths. Upper right: Scatter plot showing that there is no significant correlation between the timing residuals and 
the transit depths. Lower right: the distribution of the 60 well-measured transit depths. This distribution is biased towards deeper transits, as transits shallower 
than 0.2% were not analyzed due to the low S/N. 


sample is added. The simulated datasets are then processed 
with the same pipeline used as was used in this work for the 
K2 data. 

The main conclusions from these simulations are; (i) 
The formal uncertainties in the ‘measured’ transit depths 
(^0.05%) are indeed underestimates of the actual uncertain¬ 
ties. The simulations using a constant depth have recovered 
depths with a scatter of 0.1%. This is likely the result of a 
combination of systematic effects induced by the short dura¬ 
tion of the transits and the white noise terms. The simulations 
show that an actual scatter in the depths greater than 0.15% 
is easily recoverable. The scatter of our K2 depths is 0.2%, 
so we conclude that it is real. The ground based observations 
confirm this, (ii) The formal uncertainties also underestimate 
the true uncertainties of the transit times (defined as the scat¬ 
ter of measured times after removing a best linear trend). This 
effect is worse in the presence of starspots, but it does not de¬ 
pend on the transit depth variations. The effect can be large 
enough to explain all the scatter observed in the K2 timing 
analysis and therefore this scatter is not significant. We have 
added a systematic uncertainty in quadrature to the timings of 
3 minutes, chosen to provide a best-fit reduced of 1. (iii) 
Our treatment of stellar spots does not induce detectable depth 
variations (see Kawahara et al. 2013; Croll et al. 2015), as ex¬ 
pected, mostly because we do not have the precision to detect 
them. 


3.2. Deviation from a standard transit profile 

A fold of the activity corrected data (see section IZ3l l about 
the period we determined of 9.145704 hours is presented in 
Fig. |6| The overall crudely triangular shaped transit profile 
is the result of a convolution of the intrinsic shape and the 
LC sampling time (see also Sects. |4|and|7]i. The depth of the 
folded profile is 0.6%, and of course, this represents an aver¬ 
age of the highly variable depths. The red curve is the mean 
out-of-transit normalized flux (averaged for orbital phases be¬ 
tween -0.5 to -0.25, and from 0.25 to 0.5). Note the clear 
positive “bump” in flux just after the transit egress and the 
smaller, but still marginally significant, bump just prior to 
ingress. These are significant at the 6-cr and 2-cr confidence 
limits, respectively. 

These “bumps”, which are not normal features of exoplanet 
transits, will be important for understanding the basic nature 
of the transits (see Sect. 0. We call these features the “pre¬ 
ingress bump” and the “post-egress bump”. Even in the cases 
of KIC 1255b and KOI 2700b, the other two exoplanets which 
appear to have dusty tails, the main distinguishing feature of 
their transit profiles is a post-transit depression (Rappaport 
et al. 2012; Rappaport et al. 2014) which is attributed to a 
trailing comet-like dust tail. In addition, KIC 1255b exhibits 
a pre-ingress “bump” which has been attributed to forward 
scattering in the dust tail near the head of the dust cloud (see, 
e.g., Rappaport et al. 2012; Brogi et al. 2012; Budaj 2013; 
van Werkhoven et al. 2014). By contrast, the transit profile of 
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K2-22b has no post-egress depression, and the most promi¬ 
nent bump comes after the egress, rather than before. These 
features will be crucial to the interpretation of the dust “tail” 
in this system. 



Orbital phase 

Fig. 6. — Folded light curve of K2-22b for an orbital period of 9.14570 
hours. The folded data have been averaged into 3.7-minute bins. The em¬ 
pirical out of transit rms vaiiations in the flux are 195 ppm. As explained in 
the text there is a positive “bump” in flux just after the transit egress and a 
smaller, but still marginally significant, bump just prior to ingress. The red 
horizontal lines represent the mean out-of-transit normalized flux. The black 
horizontal bar indicates the LC time of 29.4 minutes, and gives an indication 
of the inherent temporal resolution of the light curve 


Finally, even though we believe the transit is due to a dust 
tail, as a baseline reference model we attempted to fit a stan¬ 
dard transit profile of a solid planet over a limb-darkened star 
(Mandel & Agol 2002) to our folded light curve. The high 
distortion of the transit light curve due to the 30 minute sam¬ 
pling precludes obtaining precise transit parameters, but we 
were able to constrain the scaled semi-major axis from the 
transit itself to be d/Rf, = 4.2j;Q g^, with the uncertainties es¬ 
timated from an MCMC anal ysis. In fact, from the stellar 
properties obetained in section ls!2l and Kepler’s 3rd law, we 
can make a better direct estimate of d/R^, = 3.3 ± 0.2, which 
is compatible with the transit fit. We ran a final transit model 
with a Gaussian prior on d/R^, based on the inferred stellar 
density, with mean value of 3.3 and a standard deviation of 
0.2. This, in turn, allowed us to estimate a mean K2 total 
transit duration (first to fourth contact) of 46 ± 1 minutes, an 


impact parameter of b = 0.68 ±0.06 (see TablelHi and a mean 
depth of (R'p/Rff = 0.55%, where R'p can be understood as 
the mean effective radius of the dust grains. 

3.3. Verification With Independently Processed K2 Data 

We have also used the Vanderburg & Johnson (2014) flux 
time series data set for K2-22 to check against the results of 
our own pipeline. With the application of a simple high-pass 
Alter, this independently processed light curve looks nearly 
identical to the one shown in the upper panel of Fig. [T] A 
fold of the data about the period we have determined yields 
a transit profile that is essentially the same as that shown in 
Fig. |6] including the appearance of a convincing post-egress 
bump, and a somewhat less significant pre-eclipse bump. To 
the extent that our pipeline and that of Vanderburg & John¬ 
son (2014) are independent, this is a satisfying test that our 
processing has introduced no artifacts into the lightcurve. 

4. GROUND-BASED TRANSIT OBSERVATIONS 

A total of 12 transits were obtained from the ground using 
1-m class telescopes, and the observations are summarized 
in Table [T] The first observations were taken with the 1.2- 
meter telescope at the Fred Lawrence Whipple Observatory 
(FLWO) on Mt. Hopkins (AZ) using the KeplerCam instru¬ 
ment, which has a single 4K x 4K Fairchild 486 CCD with 
a 23'. 1 X 23'. 1 held of view. We successfully observed three 
different transits with photometry in the Sloan i’ Alter. 

We also conducted a follow-up transit observation with 
the Okayama 1.88-m telescope using the NIR camera ISLE 
(Yanagisawa et al. 2006), and adopting a similar observing 
scheme to that described in Fukui et al. 2014, we studied the 
transit of the target in T-band. We only slightly defocused the 
stellar images, but due to the faintness of the target and refer¬ 
ence star, the raw counts were well within the detector’s linear 
range. The typical FWHM of the target star’s PSF was ^ 14 
pixels, which corresponds to ^ 3.5". 

We obtained two additional transits with the IAC-80 (80- 
cm) telescope at the Observatorio de IzaAsa, in the Canary 
Islands. We used the wide held CAMELOT camera, with a 
EOV of 10.4' X 10.4', with observations taken in the i-band. 
No defocussing was applied. Weather conditions were clear 
and stable through the two nights. 

A total of six transits were observed with the 0.6-m TRAP- 
PIST robotic telescope (TRAnsiting Planets and Planetesl- 
mals Small Telescope), located at ESO La Silla Observatory 
(Chile). TRAPPIST is equipped with a thermoelectrically- 
cooled 2K x 2K CCD, which has a pixel scale of 0.65" that 
translates into a 22' x 22' held of view. Eor details of TRAP¬ 
PIST, see Gillon et al. (2011) and Jehin et al. (2011). The ob¬ 
servations were obtained through a blue-blocking fllteS that 
has a transmittance of > 90% from 500 nm to beyond 1000 
nm, without any additional defocus due to the faintness of the 
target star. The procedures for the observation and data reduc¬ 
tion are similar to those described by Gillon et al. (2013) and 
we refer to this paper for further details. 

Einally, additional observations were taken in queue mode 
with OSIRIS @GTC on three different nights (see Table [1]), 
covering three complete transits of K2-22b. A spectroscopic 
time series was taken in staring mode, starting w 1 h before 
the ingress, and finishing ss 1 h after the egress. The observing 
logs are summarized in Table [T] Even though the GTC data 

http://www.astrodon.com/products/filters/exoplanet/ 
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were taken in low-resolution spectral mode, for the purpose 
of transit timing, we used the data integrated over wavelength 
(i.e., an effective white light transit). The spectral dependence 
of the transits is discussed in detail in Sect. [T] 

In all, 15 transits were measured from the ground. A log of 
these 15 observations is given in Table[T] 

Each of the above instruments has a different procedure for 
reducing the light curves, and we refer to the corresponding 
literature for a more detailed explanation (Okayama; Fukui et 
al. 2011; TRAPPIST; Gillon et al. 2013; FLWO; Holman et 
al. 2006; 1AC80: Fazaro et al. 2015; GTC: see Sec.Hli. All 
transit light curves were obtained by comparing the fluxes of 
the host star to a reference light curve made by combining up 
to several comparison stars. They were all adjusted for dif¬ 
ferential airmass corrections and a second-order polynomial 
was fitted to the out-of-transit part of the light curves to re¬ 
move long-term trends. The times of observation for all light 
curves were transformed into BID (TDB format, see Eastman 
et al. 2010) to compare them with the K2 observations. 

The 15 transit profiles measured from the ground are pre¬ 
sented in Fig. |7] The profiles were fit with the same simple 
three-segment transit model used in the K2 data analysis. The 
parameters were transit depth, time of transit, transit dura¬ 
tion, and ingress time. We also added two parameters to fit 
for any fiducial linear trends with time. In several cases a 
transit fit barely represents an improvement over a straight- 
line fit, whereas in some cases a deep transit is detected, con¬ 
firming the depth variations (see Fig. |7]i. We first scaled the 
uncertainties in the flux measurements to be equal to the stan¬ 
dard deviation of the flux residuals with respect to the best 
fit model. Correlated noise was taken into account using the 
time-averaging method, in which the ratio of the standard de¬ 
viation of the time-averaged residuals and the standard devia¬ 
tion expected assuming white noise is calculated over a range 
of timescales (Pont et al. 2006; Winn et al. 2008). In our 
case, we obtained the final values of this /3 parameter as the 
mean of the ratios with timescales from 10 to 30 minutes (see 
Table [U, and multiplied the initial uncertainty by /3 when es¬ 
timating the standard function. In those 9 cases where the 
transit fit improves the minimum by at least 35, an MCMC 
routine was used to estimate the uncertainties in the transit pa¬ 
rameters. The 9 new transit times obtained from the ground 
based observations are summarized in Table |2l We combined 
the K2 transit times with the new 9 transit times to obtain the 
final orbital ephemeris. The new transit times also show a 
high level of scatter (x^ = 30 for 9 O-C determinations). 

The O-C values for all the well measured K2 and ground- 
based transits are summarized in Fig. 0 (left panel). We used 
these data to determine a best fit orbital period and its uncer¬ 
tainty (after multiplying by the square root of the reduced x^)- 
We repeated the process of fitting the O-C points, but this time 
with a quadratic function, to set upper bounds on the deriva¬ 
tive of the orbital period. The period and period derivative 
results are summarized in Table |4] 

The extra source of scatter in the ground-based transit tim¬ 
ings could be caused by changes in the shape of the transit 
light curve (see Croll et al. 2015). In principle these tim¬ 
ing variations should be accompanied by transit duration or 
shape variations, but our light curves are not precise enough 
to allow the detection of such correlated variations (see right 
panel of Fig. |8ll. We cannot discard the possibility that part 
of the scatter is due to the use of a symmetric transit profile 
to compute transit times, when indeed some of the transits 


appear slightly asymmetric. From this set of 9 high-quality 
ground-based transit measurements we obtained a weighted 
total transit duration average of 50 ± 2 min, slightly longer 
than the duration obtained from the K2 photometry. This dif¬ 
ference of 4 ± 2 minutes, can be explained by the non-zero ca¬ 
dence of the ground-based observations, which has not been 
taken into account in these fits. A final value of 48 ± 3 minutes 
is quoted for the total transit duration in Table |4] 


5. PROPERTIES OF THE HOST STAR 
5.1. Imaging and color information 

An early inspection of the SDSS images of the target star, 
K2-22, showed that it is quite cool and likely an M star. The 
image also indicated the presence of a faint companion at a 
distance ^2" in the South-West direction. 

We observed the target star and its companion with Hyper- 
Suprime Cam (HSC, Miyazaki et al. 2012) on the Subaru 8.2- 
m telescope on 25 January, 2015 (UT). The sky condition on 
that night was clear and photometric. The HSC is equipped 
with 116 fully-depleted-type 2048 x 4096 CCDs with a pixel 
scale of 0.17". We took images of 3 s and 60 s exposures 
through g-, i-, and z-band Alters (see Fig. |9]l. The host star 
is saturated in i- and z-band with 60-s exposures, while the 
companion star is not clearly seen in the g-band image with 3- 
s exposure. We thus use the 60-s exposure images for g-band 
and 3-s exposure images for i- and z-band for the analysis 
presented here. 

The presence of a stellar companion can be problematic, 
particularly when the distance between the stars is smaller 
than the pixel size of Kepler. We measure the flux ratios of 
the companion to the host star in the HSC g-, i-, and z-band 
images (see Fig.|9]l by using the DoPHOT program (Schechter 
et al. 1993), which performs PSF-fltting photometry while de¬ 
blending stars. The background star is situated at a distance of 
1.91", with a position angle of 134 degrees west from north 
(see Table [3] for the coordinates). We obtain flux ratios of 
4.0 ±0.5%, 8.6 ±0.5%, and 11.7 ±0.5% for the g, i, and 
z bands, respectively, where the uncertainties have been in¬ 
creased to take into account systematic sources of noise. 

The contribution of the background star points toward a 
cooler and fainter companion. The bright star has been char¬ 
acterized as part of the K2-TESS catalog, in which effective 
temperatures are obtained from the colors of the stars (Stassun 
et al. 2014). The temperature reported there for the brighter 
target star is ^ 3700 K, whi ch m eans that both stars are likely 
to be dwarfs stars (see Sect. 15.2b . The colors used in that anal¬ 
ysis could have been blended due to the proximity of the com¬ 
panion. We also obtained the ugriz SDSS magnitudes and the 
J 2MASS magnitude for both stars, confirming that the col¬ 
ors used in Stassun et al. (2014) were correct. We re-derived 
the flux ratios in g, i and z band (5.4%, 6.5% and 12% re¬ 
spectively) confirming the measurements obtained using the 
Subaru images. The minimum uncertainties used in the fits 
are 0.05 magnitudes to take into account systematic effects. 
Adopting the full line-of-sight extinction of Ay = 0.17 based 
on the expected distance to the objects, and a l og g of 5 due 
to their red colors and dwarf nature (see Sect. 15.2b . we ob¬ 
tained temperatures of 3800 ± 150 K for the bright star and 
3150 ± 200 K for its fainter companion. 

The magnitudes and colors of the two stars based on a com¬ 
bination of the SDSS, Subaru, and 2-MASS photometry are 
summarized in Table |3] 
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TABLE 1 

Ground based transit observations of K2-22b 


Date 

[UT] 

Start 

End 

Epoch 

Telescope 

Instrument 
or filter 

Number of 
data points 

Median time 
between points [min] 

Aimiass^ 

ab 

[ppm] 

/Jb 

2015 Jan 15 

07:30 

12:45 

595 

1.2mFLWO 

Sloan i’ 

136 

2.3 

1.89-J- I.I5-> 1.25 

2800 

1 

2015 Jan 17 

14:40 

18:23 

601 

Okayama 

J-band 

102 

2.1 

2.32I.I8-:> 1.69 

3000 

I 

2015 Jan 23 

07:13 

12:06 

616 

1.2mFLWO 

Sloan i’ 

126 

2.3 

1.78 1.15 1.24 

3500 

I 

2015 Jan 27 

02:33 

05:55 

626 

IAC-80 

Sloan i 

II7 

1.7 

I.22-:> I.I3 1.28 

4070 

1.24 

2015 Jan 29 

03:15 

06:45 

634 

GTC 

RIOOOR 

44 

4.6 

I.I2-)- I.II -j- 1.53 

570 

1.13 

2015 Feb 02 

03:24 

09:35 

642 

TRAPPIST 

Blue blocking 

314 

1.1 

2.34-j- I.I8 1.42 

5500 

1.15 

2015 Feb 04 

02:16 

07:06 

647 

IAC-80 

Sloan i 

260 

1.1 

1.18 I.II 1.79 

4300 

1.00 

2015 Feb 12 

02:32 

06:28 

668 

TRAPPIST 

Blue blocking 

207 

1.1 

2.60-!-1.18 1.18 

4070 

1.52 

2015 Feb 14 

00:10 

03:30 

673 

GTC 

RIOOOR 

40 

5.4 

1.45 I.II 1.12 

920 

1.46 

2015 Feb 15 

03:24 

06:42 

676 

GTC 

RIOOOR 

38 

5.4 

I.I2-)- I.I2-)-2.07 

690 

1.26 

2015 Feb 18 

05:22 

12:30 

684 

I.2mFLWO 

Sloan i’ 

188 

2.3 

1.88 1.15 1.90 

2270 

1.50 

2015 Feb 23 

03:27 

07:39 

697 

TRAPPIST 

Blue blocking 

213 

1.1 

1.50-:-I.I8-:-1.30 

2600 

1.38 

2015 Feb 25 

02:27 

05:24 

702 

TRAPPIST 

Blue blocking 

137 

1.1 

I.7I 1.18 -)• 1.18 

3230 

I 

2015 Feb 26 

04:38 

08:48 

705 

TRAPPIST 

Blue blocking 

202 

1.1 

1.23 -)• 1.18 -)• 1.69 

2760 

I 

2015 Mar 21 

02:00 

05:37 

765 

TRAPPIST 

Blue blocking 

182 

1.1 

1.43 1.18 1.26 

2800 

I 


^ The airmass range is shown as zq ^min where zo and represent the airmass at the beginning and at the end of the night, respectively, and Zniin 

represents the minimum airmass. 

^ cr refers to the flux scatter respect to the best fit model, whereas ^ represents the level of correlated noise (see section[4j. 


TABLE 2 

Ground based transit times 


TABLE 4 

Properties of K2-22b 


Epoch 

Time of transit (BID) 

Uncertainty (days) 

595 

2457037.8608 

0.0007 

626 

2457049.6705 

0.0022 

634 

2457052.7269 

0.0008 

647 

2457057.6831 

0.0021 

673 

2457067.5839 

0.0015 

676 

2457068.7295 

0.0005 

684 

2457071.7778 

0.0018 

697 

2457076.7327 

0.0012 

705 

2457079.7805 

0.0012 


See section[4|for details. 


TABLE 3 

Properties of the host star and companion 


Parameter (units) 

Host star 

Companion 

RA (J2000) 

11:17:55.856 

11:17:55.763 

DEC (J2000) 

-1-02:37:06.79 

-1-02:37:05.48 

u-mag (SDSS) 

19.07 ±0.05 

21.68±0.17 

g-mag (SDSS) 

16.44 ±0.05 

19.61 ±0.05 

r-mag (SDSS) 

15.01 ±0.05 

18.79 ±0.05 

i-mag (SDSS) 

14.38 ±0.05 

17.34 ±0.05 

z-mag (SDSS) 

14.05 ±0.05 

16.34 ±0.05 

J (2MASS) 

12.74 ±0.05 

14.87 ±0.05 

H (2MASS) 

12.09 ±0.05 

14.27 ±0.05 

Ks (2MASS) 

11.91 ±0.05 

13.93 ±0.05 

Tcff (K) 

3830 ±100 

3290 ±120 

logg 

4.65±0.12 


[Fe/H] 

0.03 ±0.08 

0.06 ±0.20 

M* (Rq) 

0.60 ±0.07 

0.27 ±0.05 

R* (Mq) 

0.57 ±0.06 

0.30 ±0.08 

L, {Lq) 


u.uiu_ooo5 

Spec Type 

MOV ±1 

M4V±1 

Distance (pc) 

225 ±50 

225 ±50 


The magnitudes and colors ai'e taken from a combination 
of the SDSS and 2-MASS photometry. The stellar parame¬ 
ters are inferred from the combined analyses of the Keck- 
HIRES, IRTF-SpeX, and UH88-SNIFS spectra. See sec- 
tion[5]for details. 


Parameter 

Value 

Orbital Period^ (days) 

0.381078 ± 0.000001 

Orbital Period (hr) 

9.145872 ±0.000024 

Transit center^ (BJD) 

2456811.1208±0.0006 

Porb/Forb (yr-‘)“ 

< 3.5 X lO-’' 

Total transit duration^ (min) 

48±3 


3.3±0.2 

Impact pai'ameter, F 

0.68 ±0.06 

d(AU) 

0.0088 ±0.0008 

RpiRmf 

<2.5 ±0.4 

Mp(Mjf 

< 1.4 

ej 

17° ±1.5° 

Mdiist (g s“')® 

10 " 

£2 (leading tail only; units of R*)^ 

0.19-0.48 

^grain 

0.3-0.5 

Impact parameter, 

0.42-0.78 

Pi/p2 (two-tail model, £[ =£ 2 )^ 

<0.5 


^ Derived from the K2 and ground-based observations. 

^ Based on the mass and radius of the host star given in Table[3] 
and Kepler’s 3*^^ law. 

^ Derived from the K2 observations using a standai'd Mandel & 
Agol (2002) fit to a hard-body transit. 

^ Based on the shallowest K2 transits. 

^ 2-cr limit based on the Keck RVs. 

^ Estimate of the half angle subtended by the star at the position 
of the planet. 

® Following the type of estimate made in Appendix D of Rap- 
paport et al. (2014) 

^ 90% confidence limits based on the models of section 0 The 
leading and trailing tails are assumed to have exponential scale 
lengths and maximum optical thicknesses, £ 2 , £\, p 2 , and pi, 
respectively. 


5.2. Spectral Studies 
5.2.1. NOT-FIES spectrum 

An exploratory spectrum was obtained on 13 Feb 2015 
with the Fibre-fed Echelle Spectrograph (Frandsen & Lind- 
berg 1999, Telting et al. 2014) mounted at the 2.56-m Nordic 
Optical Telescope of Roque de los Muchachos Observatory 
(La Palma, Spain). We used the 1.3'^ Med-Res fibre which 
provides a resolving power of R = 46000 over the spectral 
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Fig. 7.— Fifteen transits of K2-22b observed from the ground, plotted to the same vertical scale. The transit models are each vertically offset by 0.02 for 
clai'ity in presentation. When necessary, the observations are binned to have a cadence close to 5 minutes. In spite of the weaker statistics for some of these, it is 
appai'ent that the transit depths vaiy considerably. 



50 100 150 200 250 40 50 60 70 

BJD -2456800 Total transit duration [min] 


Fig. 8.— Left panel: Timing residuals after removing the best fit linear orbital ephemeris. The uncertainties on the K2 timings have been increased to take into 
account systematic effects. Right panel: The uncertainties in the transit durations are too high to allow for a detection of a possible correlation between transit 
duration and transit time. 
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Fig. 9.— SubaiTj/HSC+z image of K2-22 shows a secondaiy companion 
star ~ 2" away, with a position angle of 134 degrees west from north. In this 
hand, the flux ratio between the faint and bright star is 12%. 

range 3640-7360 A. The TIES data revealed a single-lined 
spectrum. Although the low signal-to-noise ratio does not al¬ 
low us to perform a reliable spectral analysis, a comparison 
with a grid of stellar templates from (Valdes et al. 2004) and 
(Bochanski et al. 2007) confirmed that the host is a cold dwarf 
star. 

5 . 2 . 2 . Keck-HIRES spectra 

We also acquired five spectra of the host star with the Keck 
Telescope and HIRES spectrometer using the standard setup 
of the California Planet Search (CPS, Howard et al. 2010). 
Over the course of four nights, (5-8 Eeb 2015) we observed 
under clear skies and average seeing (~1.2"). Exposure times 
of less than 5 minutes resulted in SNRs of 5-10 per pixel. 
Extensive scattered light and sky emission were unavoidable 
for exposures on 5 Eeb 2015 due to the close proximity of the 
nearly full moon to the target star. All spectra were taken using 
the C2 decker, which is 0.87" wide by 14" long, resulting in 
a spectral resolution of R=60,000. The slit was oriented to 
minimize the amount of contamination from the companion 
star. 

We estimated the effective temperature T^g, surface gravity 
logg, iron abundance [Ee/H], and projected rotation veloc¬ 
ity vsin/* of the host star from the co-added HIRES spec¬ 
trum, which has a S/N ratio of about 20 per pixel at 6000 A. 
We used a modified version of the spectral analysis tech¬ 
nique described in Gandolfi et al. (2008), which is based on 
the use of stellar templates to simultaneously derive spectral 
type, luminosity class, and interstellar reddening from flux- 
calibrated, low-resolution spectra. We modified the code to 
fit the co-added HIRES spectrum to a grid of templates of 
M dwarfs recorded with the SOPHIE spectrograph (Bouchy 
et al. 2008). We retrieved from the SOPHIE archival] the 
high-resolution {R = 75000), high S/N ratio (>80) spectra of 
about 50 bright red dwarfs encompassing the spectral range 
K5-M3V. The stars were selected from the compilation of 
(Lepine et al. 2013). We downloaded spectra with no simulta¬ 
neous thorium-argon observations, to avoid potential contam- 

Available at http: //atlas . obs-hp. f r/sophie/ 


ination from the calibration lamp. 

The photospheric parameters of the template stars were ho¬ 
mogeneously derived from the SOPHIE spectra using the pro¬ 
cedure described in Maldonado et al. (2015), which relies 
on the ratios of pseudo-equivalent widths of different spec¬ 
tral features. Unfortunately, the low S/N ratio prevented us 
from directly applying this technique to the co-added HIRES 
spectrum of the target. 

Prior to the fitting procedure, the resolution of the tem¬ 
plate spectra was somewhat degraded to match that of the 
HIRES spectrograph {R = 60000) - by convolving the SO¬ 
PHIE spectra with a Gaussian function mimicking the dif¬ 
ference between the two instrument profiles. A corrective 
radial velocity shift was estimated by cross-correlating the 
observed and template spectra. We restricted the spectral 
range over which the fit is performed to 5500-6800 A and 
masked out the regions containing telluric lines. We se¬ 
lected the 5 best fitting templates and adopted the weighted 
means of their spectroscopic parameters as the final esti¬ 
mates for the target star. We found that the target star has 
an effective temperature of 7’eff= 3780 ± 90 K, surface grav¬ 
ity of logg = 4.65 ± 0.12 (logio cms““), and iron abundance 
of [Ee/H] = 0.05 ± 0.08 dex. We also set an upper limit of 
1.5 km s“' on the projected rotation velocity v sin i* by fitting 
the profile of several clean and unblended metal lines to the 
PHOENIX model spectrum (Husser et al. 2013) with the same 
parameters as the target star. Eigure [TO] shows the co-added 
HIRES spectrum in the spectal region around the Ha line, 
along with the best fitting SOPHIE template. 

5 . 2 . 3 . IRTF-SpeX spectrum 

A near-infrared spectrum of the star was also obtained using 
the updated SpeX (uSpeX) spectrograph (Rayner et al. 2003) 
on the NASA Infrared Telescope Eacility (IRTE). SpeX ob¬ 
servations were taken using the short cross-dispersed mode 
and the 0.3" x 15" slit, which provides simultaneous cover¬ 
age from 0.7 to 2.5 pm at R ~ 2000. The slit was aligned to 
capture both the target and companion spectrum. The pair 
was nodded between two positions along the slit to subse¬ 
quently subtract the sky background. Ten spectra were taken 
following this pattern, which provided a final S/N of ~ 100 
per resolving element in the X-band for the primary, and ~ 
25 for the companion. The spectra were flat fielded, extracted, 
wavelength calibrated, and stacked using the SpeXTool pack¬ 
age (Cushing et al. 2004). An AOV-type star was observed 
immediately after the target, which was used to create a tel¬ 
luric correction using the xtellcor package (Vacca et al. 2003). 

We analyzed the SpeX spectra obtained for both the bright 
and faint star, and both show strong atomic and weak CO 
absorption, as is expected for dwarf stars. Comparison with 
dwarf and giant NIR templates from the IRTE library (Rayner 
et al. 2009) rule out the possibility of either component being 
evolved. Metallicity was derived from the SpeX data using the 
procedures from Mann et al (2013), who provide empirical re¬ 
lations between atomic features and M dwarf metallicity, cal¬ 
ibrated using wide binaries. Teff was calculated using the em¬ 
pirical calibration from Mann et al. (2013), which is based on 
stars with Teff determined from long-baseline optical interfer¬ 
ometry (Boyajian et al. 2012). This analysis yielded a metal¬ 
licity of 0.00±0.08 and a Teff of 3880±85 K for the primary; 
and a metallicity of 0.06±0.20 and a Teff of 3290±120 K for 
the companion. Both the metallicity and T^g determinations 
for the primary star are consistent with those derived from the 
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Fig. 10.— HIRES co-added spectrum of K2-22 (black line) encompassing 
the Ha line. The best fitting template spectrun is overplotted with a thick red 
line. 

analysis of the HIRES spectra. 

5.2.4. UH88 SNIPS spectra 

Finally, spectra of both K2-22 and its companion star were 
obtained with the SNIPS integral field spectrograph on the 
UH88 telescope on Manna Kea during the night of UT 31 
March 2015. The two stars were spatially resolved in the im¬ 
age cubes. SNIPS spectra cover 3200-9700 A with R Ki 1000, 
do not suffer from slit effects, and have been precisely cali¬ 
brated by extensive observations of spectrophotometric stan¬ 
dards (Lantz et al. 2004; Mann et al. 2011; Mann et al. 2013). 
The wavelength coverage and resolution of SNIPS is more 
than adequate for measuring the strength of key molecular 
bands and atomic lines as indicators of Teff and gravity for 
M dwarf stars. SNR >100 was obtained for the primary. 

The effective temperature of the target star was derived in¬ 
dependently from the SNIPS spectrum by comparing it to 
Dartmouth Stellar Evolution model predictions in a man¬ 
ner that has been calibrated to retrieve the bolometrically- 
determined temperatures of nearby stars with measured an¬ 
gular radii (Boyajian et al. 2012; Mann et al. 2013). Radius 
and mass were then derived from empirical relations based 
on an expanded set of calibrator stars (Mann et al. 2015), 
with masses obtained from the Delfosse et al. (2000) sam¬ 
ple. These yield Teff = 3780 ± 60 K, R* = 0.55 ± 0.03 R©, 
M* = 0.57 ± 0.06M©, and logg = 4.72 ± 0.05, in quite good 
agreement with the other determinations of these parameters 
discussed above. 

Adopting the mean [Fe/H] and Teff derived from SpeX and 
HIRES we derived R*, M*, and L* for the primary star, and we 
did the same with the SpeX parameters for the secondary star. 
To this end we utilized empirical Teff-pe/H]-/?*, Teff-pe/H]- 
M*, and 7’eff-[Fe/H]-L* relations derived from the Mann et 
al. (2015) sample. Accounting for errors in Teff, [Fe/H] and 
the Mann et al. (2015) relations we computed physical pa¬ 
rameters that are listed in Tabled 

We also derived independent estimates of the distance to the 
stars, based on the colors and also on the spectroscopic param¬ 
eters, reaching similar values with different methods and for 
both stars. We conclude that both stars are likely bound and 
at a distance of 225 ± 50 pc. 

Constraints on the age of the system can be derived from 
the stellar rotation period using gyrochronology (see, e.g., 
Barnes 2007). It has been suggested that the classical rela¬ 
tionships may not work for M dwarfs because of the large 


spread in rotation periods of stars in clusters (e.g., Reiners 
8c Mohanty 2012). However, the ob serve d rotation period of 
the target star, 15.3 days (see Sect. 12.3b . sits slightly above 
the main period distribution of M dwarfs in the Hyades and 
Praesepe clusters (Delorme et al. 2011), possibly indicating 
an older age. On the other hand, McQuillan et al. (2013) 
finds a bimodal distribution of the rotation periods of field 
M dwarfs, likely associated with two populations with differ¬ 
ent median ages. K2-22 belongs to the shorter period (10-25 
days), younger stellar population group. Both observations 
suggest that the age of the system could be between 1 and a 
few Gyrs. 

The overall best-determined physical properties of the two 
stars, based on the spectroscopic observations and analyses, 
includingM*, R*, L*, Teff, logg, metallicity, and distance, are 
given in Table [3] 

6. RADIAL VELOCITIES 

The five Keck HIRES spectra obtained can also be used 
to place constraints on the mass of the putative planet. The 
standard CPS pipeline is used to convert from raw spectra to 
two-dimensional spectra. Each of the three HIRES CCDs are 
independently reduced with flat fields, sky subtraction, and 
cosmic-ray removal. The pixel columns at each wavelength 
are then summed, resulting in flux as a function of wavelength 
for each pixel. Consistent wavelength solutions are insured by 
aligning a carefully chosen set of Thorium-Argon emission 
lines onto the same pixels at the beginning of each night’s 
observations. 

The systemic radial velocity of each star is measured using 
the A-band and B-band telluric line features. Using the tel¬ 
luric lines as the wavelength fiducial, the relative placement 
of the stellar absorption lines is measured, and referenced to 
stars of known radial velocity (Chubak et al. 2012). These 
radial velocity measurements are made relative the Earth’s 
barycenter and are accurate to ±0.3 km/s. No RV variabil¬ 
ity in phase with the orbit of the target is detected, and the 
radial velocities have an rms of 0.3 km/s, compatible with the 
expected uncertainties. 

We used these five Keck HIRES points to set a formal l-cr 
upper limit on the RV amplitude of the host star of 280 m/sec, 
yielding a 2-a upper limit on the planet mass of 1.4 My. 
This is not highly constraining in the context of a small rocky 
planet, but it does rule out non-planetary scenarios (assuming 
that the source of the photometric dips is the brighter star). 

7. GTC MULTICOLOR OBSERVATIONS 

7.1. Wavelength Dependent Transits 

Spectro-photometric observations for 3 complete transits 
were obtained with OSIRIS on the GTC (see also Sect. [Hi. 
The GTC instrument OSIRIS consists of two CCD detectors 
with a field of view (FOV) of 7.8' x 7.8' and a plate scale of 
0.127" per pixel. For our observations, we used the 2x2 
binning mode, a readout speed of 200 kHz with a gain of 0.95 
e-/ADU and a readout noise of 4.5 e-. We used OSIRIS in 
its long-slit spectroscopic mode, selecting the grism RIOOOR 
which covers the spectral range of 520-1040 nm with a res¬ 
olution of R = 1122 at 751 nm. The observations and results 
we present here were taken using a custom built slit of 12" 
in width, with the target and a comparison star both located 
in the slit. The use of a wider slit has the advantage of re¬ 
ducing the possible systematic effects that can be introduced 
by light losses due to changes in seeing and/or imperfect tele¬ 
scope tracking (Murgas et al. 2014). 
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During the first transit observed with the GTC (29 Jan) we 
took as a reference a very close comparison star to the east of 
the target. However, for the two last transits (13 and 14 Feb) 
we took a reference star with a similar brightness to K2-22 
and located at a distance of 2.7' from the target. The posi¬ 
tion angle of the reference star with respect to the target was 
-79.2°. The two stars were positioned equidistantly from the 
optical axis, close to the center of CCD#1, while CCD#2 was 
turned off to avoid crosstalk. 

The basic data reduction of the GTC transits was per¬ 
formed using standard procedures. The bias and flat field 
images were produced using the Image Reduction and Anal¬ 
ysis Facility (IRAlj3) and were used to correct the images 
before the extraction of the spectra. The extraction and wave¬ 
length calibration was made using a PvRAlPI script written 
for GTC@OSIRIS long-slit data. This script automated some 
of the steps to produce the spectra such as: extraction of each 
spectrum, extraction of the corresponding calibration arc, and 
wavelength calibration (using the HgAr, Xe, Ne lamps pro¬ 
vided for the observations). All spectra were aligned to the 
first spectrum of the series to correct for possible shifts in the 
pixel/wavelength solution during the observations caused by 
flexures of the instrument. Several apertures were tested dur¬ 
ing the reduction process, and the one that delivered the best 
results in terms of low scatter (measured in rms) in the points 
outside of transit for the white light curve was selected. The 
results presented here were obtained using apertures of 28, 
40, and 44 pixels in width for the three transit observations, 
respectively. Final spectra were not corrected for instrumen¬ 
tal response nor were they flux calibrated. Fig. [TT] shows the 
extracted spectrum of K2-22. 



Fig. 11.— Visual explanation of our different choices for splitting the GTC 
observations into several different broad color bands. These are centered at 
652.5 nm and 820 nm for the ‘red’ and ‘blue’ colors, and 625 nm, 735 nm, 
and 850 nm respectively, for the 3-color bands. 


The universal time of data acquisition was obtained using 
the recorded headers of the spectra indicating the opening and 
closing time of the shutter in order to compute the time of 
mid exposure. We then used the code written by Eastman et 

IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron¬ 
omy (AURA) under a cooperative agreement with the National Science Foun¬ 
dation. 

Python environment for IRAF 


TABLES 

Multicolor observations summary 


Date 

SsiO-eWnm 

<i680-790nm 

'5790-900nm 

a 

Jan 29 

0.50 ±0.06% 

0.47 ±0.05% 

0.46 ±0.05% 

0.13±0.55 

Feh 14 

0.48 ±0.04% 

0.49 ±0.03% 

0.44 ±0.03% 

0.11 ±0.37 

Feb 15 

0.95 ±0.04% 

0.83 ±0.03% 

0.72 ±0.04% 

0.83 ±0.23 


These transit depths are measured from the data alone, and do not take into 
account contamination from the faint companion or airmass connections. 
The calculated Angstrom exponent, a, is, however, corrected for the small 
dilution factor of the faint companion star. 


al. (2010010 compute the BID time using the mid-exposure 
time for each of the spectra to produce the light curves ana¬ 
lyzed here. 

We constructed light curves over several spectral ranges to 
search for color dependencies of the transit shape and depth. 
In the case of the white light curve, the flux was integrated 
over almost the entire wavelength range of the observed spec¬ 
tra, between 570 and 900 nm, but avoiding the blue and red 
ends of the spectra where the SNR is lower. This step is 
particularly important, because observations at redder wave¬ 
lengths suffer from fringin^^ whereas using the data at bluer 
wavelengths not only did not improve the quality of the light 
curves, but increased the occurrence of outliers. Three nar¬ 
rower color light curves were created from the raw spectra for 
each star, by integrating all counts over the spectral ranges 
570-680 nm, 680-790 nm, and 790-900 nm (see Fig. fllT) . 
The final light curves used in the analyses were created by 
dividing the K2-22 light curves by the comparison star light 
curves. The scatter (std) for the relative white light curves can 
be found in Table [T] 

The white light transit profiles of the three separate GTC 
transits were shown earlier in Fig. |7] It is obvious that the 
depth of the transits change from one night to the next, con¬ 
firming the rapidly varying nature of this object. Note that 
between transits 2 and 3, only 27.44 hours had passed and the 
transit depth nearly doubled. 

To further investigate the nature of these changes in the 
GTC transit curves, in Fig. [12] we plot the three color light 
curves for each transit separately. Transit 3 (bottom panel) is 
the deepest and shows a clearly increasing transit depth to¬ 
wards the blue. The SNR is lower in the other two sets of 
color light curves and the transits are shallower, and thus there 
is no clear trend in transit depth with wavelength. In order 
to interpret quantitatively the color dependence of the depths 
in the third GTC transit, we fit our standard three-segment 
transit profile model to the data from each of the nights, with 
three different transit depths to evaluate how the transit depth 
changes with color. We also fit each transit with five addi¬ 
tional parameters, to construct a polynomial that reaches sec¬ 
ond order in time from mid-transit and it linearly depends on 
the mean-subtracted airmass and full-width half maximum of 
the images of the host star. The total number of parameters 
is 21 for each observation, with a total 120, 132 and 114 data 
points in each of the three nights. After finding the best-fit 
model, we set the error of the flux measurements of each light 
curve to provide a best-fit standard value equal to the num¬ 
ber of degrees of freedom. We then took correlated noise into 
account computingthe (3 factors for each light curve and each 
night (see section |4]i, and multiplying the errors by their cor- 

http://astroutils.astronomy.ohio-state.edu/time/ 

http://www.gtc.iac.eS/instruments/osiris/osiris.php#Fringing 
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Fig. 12.— GTC observations split among three color bands. On the left panels we show the raw light curves, which have been median-normalized. The 
observations shown on the right hand side have been corrected for time, airmass and seeing effects. The first two transits were shallower, and depth variations are 
hard to notice, but during the third night, when the white light transit is deepest, color variations are clearly obseiwed. 


responding /3 factor. We use an MCMC routine to obtain pos¬ 
terior distributions for the transit depths, which are marginal¬ 
ized respect to all the other 18 model parameters (see Table|5]l. 

From the marginalized posterior distributions we compute 
a quantity called the “Angstrom exponent”, a, which is de¬ 
fined as -c/In cr/r/In A , where tr is the effective extinction 
cross section. Here we take the transit depth (on each night) 
at wavelength. A, to be proportional to the effective cross sec¬ 
tion at A (under the assumption that the dust tail is optically 
thin). In order to use the transit depths, we first multiply each 
individual transit depth by (1 -fD,) to apply a dilution correc¬ 
tion, where each D, is the flux ratio between the faint com¬ 
panion and the host star evaluated at the flux weighted center 
of each selected band. Thes e D, values are obtained from our 
SED models (see Sect. EB, and have values of 4.4%, 6.6% 
and 8.9% with increasing wavelength, with uncertainties of 
0.2%. We then find a = 0.83 ± 0.23, for the transit observed 
on February 15. For the other two transits, the values of a 
are 0.11 ±0.37 and 0.13 ±0.55, both slightly positive, but 
not statistically significantly different from zero. We made 
no corrections for the wavelength-dependent limb darkening 
properties, but utilized the quadratic limb darkening coeffi¬ 
cients of Claret & Bloemen (2011) to estimate that the val¬ 
ues of a would be lowered by between 0.15 and 0.02, for an 
equatorial transit (which is unlikely) and an impact parameter 
of b = 0.7, respectively. For higher impact parameters, up to 
0.85 (see Table |4]i, the value of a could actually be raised by 
up to 0.15, and become even more significantly different from 
zero. 

7.2. Interpretation in Terms of Dust 

We interpret the value of the Angstrom exponents in terms 
of Mie scattering (for spherical dielectric dust particles) with a 
variety of different compositions. We computed the Angstrom 
exponent, a, over the range 630-840 nm for a set of power-law 
distributions for the dust particle sizes, with dN/da oc a~^. In 
order to guarantee that the total cross section converges, we 
also need to specify a maximum grain size, amax- For the sake 
of specificity we adopted an illustrative dust composition of 
corundum, but the conclusions we draw are the same for a 
number of other common refractory materials, and in fact, 
for any material with real and imaginary indices of refrac¬ 


tion of « ~ 1.6 and 0.001 < k < 0.03. We plot the computed 
Angstrom exponent in Fig. [13] as a function of flmax for five 
different power-law exponents, F. As can be seen from the 
figure, values of a in the range of ~0- 1, as indicated in Ta¬ 
ble |5| correspond to non-steep power-law indices of F ~ 1-3 
and maximum particle sizes of ^0.4-0.7 pm. In turn, these 
correspond to “effective particle sizes” of 0.2 pm to 0.4 pm, 
where the effective radius is the average grain radius weighted 
by both the size distribution and the cross section, i.e., 

P ^max P ^max 

Oef[= / a^~^ a{a,X) da / / a~^a{a,X)da (1) 

Jo Jo 

where a{a,X) is the wavelength dependent Mie extinction 
cross section for a particle of radius a. 



Maximum radius («m) 


Fig. 13.— The Angstrom exponent computed from Mie scattering over 
the wavelength range 0.63 —0.84/rm as a function of the maximum particle 
size, for five different indices, F, of the power-law grain size distribution, i.e., 
dN/da oc a~^. We adopted an illustrative particle material of corundum. 


8. WHY THE BRIGHTER STAR IS THE SOURCE OF THE TRANSITS 

Here we summarize why we are confident that it is the 
brighter star that is the source of the observed transit s. Fi rst, 
the fainter star is redder than the bright star (see Sect. ED; in 
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particular the contribution from the faint star to the total flux 
increases by a factor of 2 from r to i band and a factor of 4 
from r to z band. If we assume that the faint star is the host, 
and assume an achromatic mean transit depth, the changes in 
flux ratio should translate into correspondingly larger transit 
depths toward the red due to the dilution effect. In section |7] 
we showed simultaneous GTC transit observations in bands 
similar to r and a combination of the i and z bands. The tran¬ 
sit depth, at least on one occasion, is actually considerably 
shallower in the latter redder band, strongly suggesting that 
the bright star is the one being transited. 

From another line of argument, we note that the transit 
depths are sometimes as large as 1.3%, and if the fainter star 
(5% of the flux in g-band) were the source of the transits, it 
would have to be attenuated by 24% of its flux. Since the 
transits are highly variable in depth, and the transit profile is 
not that due to a conventional hard-body transit, it is almost 
certainly due to a ‘soft’ attenuator such as a dust tail. But, the 
fainter star has a radius of 0.3 Rq, and therefore an occulting 
dust tail would have to either cover the entire star with an op¬ 
tical depth, T, of ~0.25, or cover ~25% of the star (i.e., a tail 
thickness of ^0.08 Rq) with t ~ 1. Simulations of the dust 
tail i n KI C 1255b (Rappaport et al. 2012) and this object (see 
Sect. I9.21 i indicate that the vertical thickness of the dust tail is 
only ^0.03 Rq. This would make the scenario of a dust tail 
covering over 1/4 of the area of the fainter star highly unlikely. 

9. THE DISINTEGRATING PLANET HYPOTHESIS 
9.1. Evidence for a Disintegrating Planet 

When taken together, all these observational results point 
clearly toward a planet in a 9-hour orbit that is disintegrat¬ 
ing via the emission of dusty effluents. The lines of evi¬ 
dence pointing in this direction include: (1) erratically and 
highly variable transit depths (see Sect. |3li; (2) transit profile 
shapes from ground-based observations that are likely vari¬ 
able (though with lesser confidence that the depth changes; 
see Sects. 0] and IT); and (3) an average transit profile from 
the K2 data that exhibits clear evidence for a post-transit 
“bum p” a nd also weaker evidence for a pre-transit “bump” 
(Sect. lT^ . The first of these is highly reminiscent of the dis¬ 
integrating planet KIC 1255b (Rappaport et al. 2012), while 
the variable transit shapes are also detected in KIC 1255b 
from ground-based studies (R. Alonso et al. private communi¬ 
cation; Bochinski et al. 2015). The transit profile of K2-22b, 
with a post-transit “bump”, is different from the transit pro¬ 
files of KIC 1255b and KOI 2700b which show a post-transit 
depression as opposed to a post-transit bump. The first two 
of the above listed features point to obscuration by dusty ef¬ 
fluents coming from a planet, while the third property needs 
to be explained in this same context. In the following sec¬ 
tions we explore the significance and the interpretation of the 
transit profile. 

9.2. Quantitative Model for the Leading Dust Tail 

A dust tail emanating from a planet, as inferred in the cases 
of KIC 1255b and KOI 2700b, trails the planet as is illustrated 
in Fig. 6 of Rappaport et al. (2012; see also the middle panel 
of our Fig.fTSll. Such dust tails are the way they would be seen 
in the reference frame of the planet, and note that the motion 
of the planet is implicitly in the opposite direction from the 
tail. In this case we would say that the tail “trails the planet”. 
The reason for this is that the radiation pressure acting on the 
dust forces it into an eccentric orbit with its periastron located 


at the point where the particle was released. This orbit has a 
larger semimajor axis than that of the planet (see Appendix B 
of Rappaport et al. 2014). In turn, the larger orbit has a lower 
orbital frequency, and the particles appear to trail behind the 
planet in a comet-like tail. 

An observer viewing this system is effectively moving 
counterclockwise. Thus, the ingress to the transit is sharp as 
the “comet head” moves onto the stellar disk (see Fig. 6 in 
Rappaport et al. 2012). The trailing tail would lead to a de¬ 
pression upon egress as the tail slowly moves off of the stellar 
disk. In that case, the pre-transit “bump” is caused by forward 
scattering by the densest regions of the dust which have not 
yet reached the stellar disk. This is the case we believe we see 
in KIC 1255b (Rappaport et al. 2012; Brogi et al. 2012; Bu- 
daj 2013; van Werkhoven et al. 2014). If the tail of the planet 
is sufficiently short (i.e., compared to the radius of the host 
star) then there would be both a pre-transit and a post-transit 
“bump”, the latter of which would dominate over the rela¬ 
tively shallow post-transit depression (see also Budaj 2013). 
In this case, the pre-transit bump would be somewhat larger 
due to the asymmetry in the direction of the tail. 

A logical first guess as to how to produce a post-transit 
“bump” on the transit curve would be to reverse the direction 
of the comet-like tail. However, as we have seen, substan¬ 
tial radiation pressure forces inevitably lead to a trailing tail. 
Then, the question becomes how to produce a “leading dust 
tail” to the planet. One way to have dusty material lead the 
planet, i.e., moving faster, would be to have it overflow its 
Roche lobe (or. Hill sphere of influence) and fall in toward 
the host star. At first consideration this doesn’t seem to work 
since a rocky planet in a 9-hour orbit will not be close to Ailing 
its Roche lobe. Rappaport et al. (2013) showed that the criti¬ 
cal density for Roche-lobe overflow is largely a function of its 
orbital period, and is nearly independent of the properties of 
the host star. In particular, Eqn. (5) in Rappaport et al. (2013) 
suggests that critical density for a planet to be Ailing its Roche 
lobe is 

/11.3hry , 3 , 3 ,,, 

Peril p ^ I g/cm~1.5g/cm (2) 

where the right-hand value is for a 9-hour planet. It seems 
very unlikely that a planet with this low a mean density would 
be disintegrating via dusty effluents. Turning the problem 
around, we can ask what fraction of the Roche-lobe radius 
would be occupied by a planet with a mean density in the 
range of 5-8 g/cc, which might be more appropriate for a dust 
emitter (see, e.g., Rappaport et al. 2013). It is then evident that 
for densities which are ^3 - 5 times higher than their critical 
densities, their radii are not even factors of ^2 times smaller 
than their Roche lobes. 

Therefore, we come to the conclusion that planets with 
rocky compositions in a 9-hour orbit are underfilling their 
Roche lobes by only a factor of ~2. This is more than suffi¬ 
cient to prevent the planet from directly overflowing its Roche 
lobe. However, the potential difference between the planet’s 
surface and the Roche lobe is about half the potential differ¬ 
ence to infinity. Thus, if the mechanism that drives off the 
dust or the heavy metal vapors that condense into dust, e.g., 
via a Parker-type wind (Rappaport et al. 2012; Perez-Becker 
&. Chiang 2013) imparts the full escape velocity of the planet 
or more, then the material can certainly reach the surface of 
its Roche lobe. 

If the material reaching the Roche surface feels a substantial 
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Fig. 14.— Calculated values of j3, the ratio of radiation pressure to grav¬ 
itational forces, as a function of particle size. The blue and green curves 
(nearly superposed) are for KIC 1255b and KOI 2700b, respectively. The 
red curves are values of 0 for EPIC 201673175, for different materials. For 
KIC 1255b and KOI 2700b we adopted « = 1.65 and k = 0.01 for the real and 
imaginary parts of the index of refraction. The solid red curve is for the same 
indices. The dot-dashed curve is for the same real index but with k = 0.02, 
while the dotted, dashed, and dot-dot-dashed curves are specifically for iron, 
corundum and forsterite, respectively. The horizontal black dotted line indi¬ 
cates the value of /3 = 0.02 below which most dust particles would go into a 
leading tail. 


radiation pressure, it will still be blown back into a comet-like 
tail. By “substantial” we mean that /3, the ratio of radiation 
pressure forces to gravitational forces, exceeds a certain value 
such as /3 > 0.05. For sufficiently small values of /3, however, 
particles initially directed toward the host star (which seems 
reasonable since the gas or dust emission should commence 
on the heated hemisphere) will largely fall toward the host star 
until the grains are pushed into orbit by coriolis forces. These 
particles will form a “leading tail” in the sense that the motion 
will be in front of the planet. 

The values of /3 for the particles are computed according to 
^ L*<J ^ (g(a,A)) 

AncGM^:fj, l^TTcGM^pa 

where L* and M* are the luminosity and mass of the host star, 
a is the particle radius, a is the particle cross section at wave¬ 
length A, {a) is the dimensionless cross section in units of 
TTo^ averaged over the stellar spectrum (see, e.g., Kimura et 
al. 2002), and p and p are the mass and material density of 
the dust grains. Since, for low-mass main-sequence stars the 
luminosity scales roughly as M^, we hnd that for a hxed ma¬ 
terial in the grains, j3 scales as 


/3 cx 


Ml (g(a,A)) 

a 


(4) 


We evaluate the dimensionless spectrum-averaged cross sec¬ 
tion with a Mie scattering code (Bohren & Huffman 1983) for 
different assumed material indices of refraction (see Croll et 
al. 2014). Plots of /3(a) for several different assumed compo¬ 
sitions are shown in Fig. [14] They are also compared against 
representative /3(a) curves for KIC 1255b and KOI 2700b. 
The latter two systems have /3(a) curves that are essentially 
a factor of 2 higher than for K2-22b, largely due to the lower 
luminosity per unit mass of the latter host star. 

Expression (|4]i raises the interesting prospect that for low 
mass stars, the value of /3 will be sufficiently small that out¬ 
flowing particles passing through their Roche lobes will be 
immune to radiation forces and act much in the same way as 



Fig. 15.— Simulated dust tails. The dust particles are launched with some¬ 
what more than the escape speed from the planet within a 30° cone centered 
on the direction of the host star. The particle sizes are chosen via Monte Carlo 
means from a power-law distribution, and the value of ^ for each particle is 
calculated from the system parameters (see text). There are 50,000 particles 
in each simulation, and each one is assumed to have an exponentially decay¬ 
ing cross section (see App. C of Rappaport et al. 2014) with a time constant of 
10"^ s. Only radiation and gravitational forces ai‘e included. The color coding 
is proportional to the logarithm of the dust particle density with white-red the 
largest to blue-puiple the lowest. The dust tails ai'e shown in the rest frame 
of the orbiting planet (implicitly moving downward in the frame). The top 
panel is a view of the orbit and dust tail from the orbital pole, and clearly 
shows both a leading and a trailing tail. The bottom panel is a view from the 
orbital plane as the dust emitting planet crosses the disk of the host star. The 
middle panel results from arbitrarily multiplying each calculated value of p 
by a factor of 4, the net effect of which is to eliminate the leading tail. 
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Fig. 16.— Schematic transit profiles for a range of dust tail parameters. 
6 is the orbital phase of the planet, and 0* is the half angle subtended by 
the host star as seen from the planet. Top panel: Trailing dust tail only. The 
exponential tail length, £, is color coded and is expressed in units of the radius 
of the host star. No forward scattering is included. Note how the shorter the 
tail length, the smaller is the post-transit depression. Middle panel: Same 
as for the top panel, but in this case a contribution from forwai'd scattering 
has been included. There is now a noticeable pre-transit bump due to the 
forward scattering, but only for the case of short dust tails is there a bai'ely 
perceptible post-transit bump; it is typically suppressed by the post-transit 
depression. Bottom panel: Here, there is both a trailing and a leading dust 
tail. The exponential scale lengths, i, of both dust tails are assumed to be the 
same, but the leading tail is taken to have twice the dust density, p, at every 
corresponding angular distance, as the trailing tail. For long tail lengths there 
are no pre- or post-transit bumps due to forward scattering, as they are pulled 
below the discernible level by the corresponding depressions caused by direct 
attenuation from the tails. By contrast, for the case of short tail lengths, i.e., 
< O.lRstair there is both a pre- and a post-transit bump, with the latter being 
slightly larger. 


Roche-lobe overflowing material. 

We next carried out a large number of simulations of dust 
particles ejected from a planet and initially moving in vari¬ 
ous directions with different speeds. For most of our simu¬ 
lations, we took the direction of the particles to be uniformly 
distributed within a cone of 30° radius and centered in the 
direction of the host star. The velocities were somewhat arbi¬ 
trarily taken to be one half the escape speed from the surface 
of the planet’s Roche lobe to infinity in the absence of the host 
star, i.e., where Mp is the planet’s mass and Ri is 

its Roche-lobe radius. This is the minimum that can be con¬ 
templated for a Parker-wind type outflow. For each dust grain, 
a radius, a, is chosen from an assumed particle size distribu¬ 
tion via Monte Carlo methods. We considered a power-law 
differential size distribution for the dust grains with a slope of 
-2, as illustrative, with the maximum and minimum sizes in 
this distribution taken to be 1 pm and 1/20 pm, respectively. 
The luminosity of the host star is taken to be that of a main- 
sequence M-K star of 0.6 Mq (see Table O. For purposes of 
the numerical calculations of the particle dynamics only, we 
used a simple analytic fit to the plots shown in Fig.[T4]of the 
form: 


P oc 


(g(a,A)) 

a 


Cl +C2a^ 
1 +€ 2 , 0 .'^ 


(5) 


where the c’s are constants to be fit, and which depend on L*, 
M*, and the dust composition. 

The results of our particle dynamics simulations are shown 
in Fig. [15] The top panel represents the trajectories of all par¬ 
ticles regardless of their size or corresponding value of /3. The 
dust density is assumed to decay exponentially in time, with a 
time constant of 10"^ sec. Note that there is a forward moving 
“leading tail” in addition to a bifurcated trailing tail. The inner 
of the two trailing tails arises from particles with small values 
of P that are launched with a substantial velocity component 
in the forward direction. In the bottom panel of Fig. [15] we 
see the two tails from the perspective of an observer viewing 
an equatorial transit; there are far more particles in the lead¬ 
ing tail. For comparison, we show in the middle panel a case 
where the computed value of P was arbitrarily multiplied by a 
factor of 4; all other parameters remained the same. Note that 
the result is a purely trailing dust tail as one would have in a 
solar system comet. This latter exercise ensures that all the 
particles have a P that exceeds a critical value of ~0.02, guar¬ 
anteeing that the particles will go into a trailing tail. The bot¬ 
tom line is that dust-emitting planets around luminous stars 
should have predominantly trailing tails while the reverse is 
true for low-luminosity stars. 

In all of these calculations, we have ignored possible ram 
pressure forces on the dust grains due to a stellar wind from 
the host star. For a justification of why this may be a good 
approximation, see Appendix A of Rappaport et al. (2014). 


9.3. Model Transits 
9.3.1. Idealized Transit Models 

The transit profiles expected for a planet with a trailing 
comet-like tail are already sufficiently complicated compared 
to a conventional hard-body planet transit. They include such 
issues as the attenuation caused by the (unknown) absorp¬ 
tion profile, typically characterized by an exponential scale 
lengtlQ possible forward scattering which is influenced by 

For a discussion of why the dust tail may fall off approximately expo¬ 
nentially, including the effects of dust sublimation on a timescale of hours, 
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green and blue curves represent the direct absorption and forward scattering 
components, respectively. Note how the forward scattering causes the peaks 
on either side of the transit (see also Budaj 2013). 

the grain size and composition (which may be changing along 
the tail as different materials sublimate at different rates); con¬ 
volution of the scattering phase function with the finite angu¬ 
lar profile of the host star; and another convolution with the 
density profile of the tail. When adding the possibility of both 
a trailing and leading tail, as we have postulated in this work, 
the situation quickly becomes even more complicated. In an 
attempt to keep the numbers of free parameters to a minimum, 
we have adopted the following seven-parameter model: an ex¬ 
ponential attenuation profile of the tail in both directions, each 
with its own scale length, l\ a relative density between the two 
tails, p \/ pi ', the size of the dust grains, a; and an overall nor¬ 
malization factor that yields the correct mean transit depth. 
In addition to these parameters associated with the dust tails, 
there is also an impact parameter for the transiting planet, and 
the time of orbital phase zero. We assume that the hard body 
of the planet itself does not make a significant contribution to 
the transit profile. 

We illustrate in Fig. [16] what some of the possible transit 
profiles might look like with just two of these parameters 
varying. First, we consider only a trailing dust tail, and we 
vary the exponential scale length (£, in units of the host star’s 
radius). These transit profiles are shown in the top panel of 
Fig dSI The tails’ exponential scale length, I, is color coded 
to help in distinguishing the different curves. No forward 
scattering is included in this case. Note how the shorter ^ is, 
the smaller is the post-transit depression. The latter becomes 
quite noticeable when ^ for the tail becomes greater than about 
30% of the radius of the star. In the middle panel of Fig. [T6| 
the contribution from forward scattering by relatively large 
particles (~0.5 pm) has been included. In all cases there is 
a noticeable pre-transit bump due to the forward scattering, 
but only for the case of short dust tails (i.e., small i) is there 
a perceptible post-transit bump; it is typically suppressed by 
the post-transit depression as the dust tail is in the process of 
moving off the stellar disk. The grain sizes assumed here are 
sufficiently large so that A/27ra, the characteristic Mie scat¬ 
tering angle, is comparable with the angular size of the host 

see Appendix C of Rappaport et al. (2014). 


star as seen from the planet. In that case, the effective angular 
scattering pattern from a single grain is somewhat larger than 
that of the angular size of the host star (i.e., ^17°). In the bot¬ 
tom panel of Fig.[T6|the transit profiles are for the case where 
there is both a trailing and a leading dust tail. The exponen¬ 
tial lengths of both dust tails have been taken to be the same 
(i.e., f 1 = ^ 2 ), but the leading tail is taken to have twice the op¬ 
tical thickness everywhere as the trailing tail. For large i there 
are no pre- or post-transit bumps, as they are pulled below the 
discernible level by the corresponding depressions caused by 
the tails. By contrast, for the case of short tail lengths, i.e., 
i < 0.2f?,, there is both a pre- and a post-transit bump, with 
the latter being slightly larger. 

9.3.2. Model Fit to the K 2 - 22 b Transit 

Finally, we have attempted to fit a simple two-tail dust 
model to the observed transit profile of K2-22b. Once the 
model and its free parameters were selected, we computed the 
scattering pattern for the dust particles (taken to be of a sin¬ 
gle fixed size throughout both tails) via a Mie scattering code 
(Bohren & Huffman 1983). Then the scattering pattern was 
convolved numerically via a 2D integral with the radiation 
profile coming from the host star (including limb-darkening). 
In turn, this net effective scattering profile was convolved in 
ID with the assumed exponential tails to produce the scatter¬ 
ing pattern as a function of orbital phase. The attenuation of 
the beam from the host star was simply computed by placing 
a double-sided exponential profile in front of a limb-darkened 
stellar disk in different longitudinal locations (i.e., as a func¬ 
tion of orbital phase) and at different vertical locations to get 
a best estimate of the impact parameter b (in units of the stel¬ 
lar radius). The sum of the absorption and forward scattering 
contributions was then added and convolved with the Kepler 
long cadence integration time. 

We first explored a wide range of parameter space, espe¬ 
cially focusing on the physically interesting parameters: b, 
the impact parameter, i\, £ 2 , pijpi, the exponential scale 
lengths and ratios of optical thickness of the two tails, respec¬ 
tively, and a, the grain size. We also considered models with 
a power-law distribution of grain sizes, but these produced 
no better fits than using a single grain size. To simplify the 
search, we fixed £i = £2 under the assumption that the dust 
sublimation time and the speed away from the planet is sim¬ 
ilar in the two tails. However, we allowed pijP 2 to be a free 
parameter since the amount of dust in the two tails could be 
quite different. From this broad search, we were able to draw 
several interesting conclusions. (1) The best fits, by only a 
small margin, are found for the single (leading) tail models, 
i.e., no trailing tail is needed. (2) For single-leading tail mod¬ 
els, the impact parameter is constrained to be 0.42 < b < 0.78 
and 0.19 < £2 < 0.48/?* (both 90% confidence limits). (3) For 
two-sided tail models the allowed ratio of optical thicknesses 
is constrained to be py/p 2 < 0.5 (assuming £\ = £ 2 ). (4) The 
best-fitting single grain-size models have 0.3 < a < 0.5 pm. 
The final parameter search and error estimation were done 
with an MCMC routine. These results are summarized in Ta- 
blelD 

The overall best-fitting model is for b = 0.65 and £2 = 
0.32/?*. The results are shown in Fig. (TT] The black his¬ 
togram is the same average transit profile as is shown in Fig.|6| 
The red curve is the overall fit to the transit profile, includ¬ 
ing the convolution with the LC integration time. The overall 
model transit profile is comprised of two components: the di¬ 
rect attenuation of the flux from the host star (green curve) 
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and the forward scattering from the dust tail back into the 
beam directed toward the observer (blue curve). Note how 
the forward scattering produces small bumps both before and 
after the transits, with the larger bump following the egress. 

9.4. Inferred Dust Sublimation Timescale 

The dust sublimation timescale depends on the equilibrium 
temperature (Teq), the mineral composition, and the size of 
the dust grains (see, e.g., Kimura et al. 2002; Appendix C 
of Rappaport et al. 2014). Teq for the dust grains in K2-22b 
is necessarily somewhat uncertain, ranging between ~1500 
K and 2100 K, depending on exactly how it is computed (to 
within unknown factors of order unity). However, if we take 
Teq = 1700 K as a typical grain equilibrium temperature, then 
the lifetime of a 1 micron grain composed of Fe, SiO, or 
fayalite is a matter of a few tens of seconds. On the other 
hand, minerals such as enstatite, forsterite, quartz, and corun¬ 
dum, would have sublimation lifetimes of 1/2, 3, 15, and 30 
hours, respectively. SiC and graphite could last for a very 
much longer time than any of these. (All the dust parame¬ 
ters for the above estimates, including the heat of sublima¬ 
tion, are taken from the conveniently tabulated values of van 
Lieshout et al. 2014). Thus, there is no one “expected” dust 
sublimation timescale that we can anticipate. If on the other 
hand, we believe that the dust tail is only a few degrees of 
the orbit in length (see Table |4]i, and that (3 k, 0.05, then the 
“inferred” lifetime is ^ 1 hour (see Eqn. (4) and (6) of Rap¬ 
paport et al. 2014 for an explanation), and therefore the dust 
composition might resemble enstatite or forserite. However, 
one should not draw too many conclusions from this since 
the actual tail length depends on a combination of numerous 
(uncertain) parameters such as the particle sizes (and indeed 
the size distribution), the parameter /3, the actual equilibrium 
temperature, and so forth. Nonetheless, we can say that min¬ 
erals such as enstatite or forsterite do seem consistent with the 
observations. 


9.5. Dust-on-Dust Collisions 

We note that there is an interesting possibility for the dust 
streaming from the planet to interact with dust that has already 
been orbiting for a substantial while - provided that the subli¬ 
mation lifetime will allow the dust to survive for a sufficiently 
long time. The time required for orbiting dust in a trailing 
dust tail (i.e., ‘outer-track orbits’ with 0.05 ^ fi > 0.02) to 
meet up again with the planet is approximately Porb/(2/3) (see 
Eqn. (4) of Rappaport et al. 2014); for K2-22b this amounts to 
10-25Porb (or r^4- 10 days). Eor particles with (3 < 0.02 the 
orbits likely lie inside the planet’s orbit (‘inner-track orbits’) 
with radii between ^94% and 98% of d. This leads to synodic 
orbital periods for those dust particles of between ^10 and 30 
Porb, fairly similar to the outer-track orbits. Thus, any orbit¬ 
ing dust that might collide with newly emitted dust must last 
for 4- 10 days. The dust grains on the inner-track orbits have 
speeds in the rest frame of the planet that range from 0-6 km 
s“' (prograde), while the corresponding range of outer-track 
orbits is 0-25 km s“' (retrograde). Thus, the inner track par¬ 
ticles would catch up to, and collide with, newly emitted dust 
grains with only a few km s“* relative velocity. By contrast, 
the outer track particles could collide with newly emitted dust 
at relatively high speeds of more than 10 km s“'. At such 
speeds, a collision would deposit a mean energy of ^1 eV per 
atomic mass unit within the dust grain. This seems likely to 
completely destroy most dust particles that happen to collide 


on the outer track. 

Additionally, any orbiting dust grains which are moving 
quickly may catch up, and collide, with other more slowly 
orbiting dust. Again, the relative speeds of any colliding dust 
particles on the inner track would likely be of order a km s“' 
and would not necessarily destroy the dust, whereas the re¬ 
verse is true for the outer-track dust grains. 

We can make a crude estimate of the collision frequency 
of a dust grain that is orbiting the host star. Take the mean 
grain number density along the observer’s line of sight during 
a typical transit to be no, and the collision cross section to be 
tr. In that case, the collision mean free path is £ = l/{noa). 
If we take the radial thickness of the dust tail that is causing 
the transits to be Ar, then Ar ~ 1 / (noa) if both the optical 
cross section for extinction and the collision cross section are 
approximately equal (i.e., the geometric cross section appro¬ 
priate for larger particles), and if the optical depth of the dust 
cloud during transits is of order unity. (The latter is necessary 
since the dust tail is thin in the vertical direction and covers 
only a small fraction of the disk of the host star (see middle 
panel of Eig. fTSl l. Using the above expression, we can esti¬ 
mate the mean collision time, Tcoh, for particles orbiting in the 
dust disk 


"^coll 


1 Ar Vorb no 
lit d Vrei (n) 


( 6 ) 


where v'orb and v^ei are the mean orbital speed in inertial space 
and the mean relative speeds of the particles in a dust disk, and 
(h) is the mean density in dust grains around the azimuth of 
the dust disk. These factors are extremely uncertain, but as an 
illustrative example, if we take iSr/d sa 0.05, v'orb/t'rei ~ 10, 
and no/(n) sa 100, then Tcoii is of order several times Porb- 

What is the fate of dust particles that collide but are not 
destroyed in the collision? There could be locations in the 
dust disk where the dust is slowed down by the collisions and 
the density builds up. In principle, these regions of enhanced 
density could be approximately fixed in the rotating frame of 
the planet and host star. However, since we see only a single 
transit-like feature during an orbital period, we conclude that 
such pile-up of dust is not significant anywhere, with the pos¬ 
sible exception of near the planet where the dust density is the 
highest. 

A careful examination of what the effects are of dust-dust 
collisions in this system or, for that matter, KIC 1255b and 
KOI 2700b are much beyond the scope of this paper, and we 
leave that study to another work. 


10. SUMMARY AND CONCLUSIONS 

We have reported the discovery in the K2 Eield-1 data of 
a new planet that is likely the third example of a planet dis¬ 
integrating via the emission of dusty effluents. The planet is 
in an ultrashort 9.1457 hour orbit about an M star. The ev¬ 
idence we presented for the presence of a dust tail includes 
erratically and highly variable transit depths ranging from 
< 0.14% to 1.3% in both the K2 data as well as in follow-up 
ground-based observations from five different observatories. 
The folded orbital light curve from the K2 data exhibits a clear 
post-egress ‘bump’ and a much less significant, but plausible, 
pre-ingress ‘bump’, that are not found in conventional hard- 
body transits. There is no post-egress depression in the flux 
as was seen in KIC 1255b or in KOI 2700b. On at least one 
observation with the GTC, the transit depths were distinctly 
wavelength dependent with the transits ~25% shallower at 
840 nm than at 630 nm. While this goes in the same direc- 
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tion as limb-darkening effects (see, e.g., Knutson et al. 2007; 
Claret & Bloemen 2011), we argue that the magnitude of the 
limb-darkening effect explains only a small fraction of what 
is observed. Furthermore, the variable behavior of the color- 
dependent transit depths strongly suggests an origin in dust 
scattering. This requires relatively non-steep power-law par¬ 
ticle size distributions with F ~ 1 to 3 with maximum sizes in 
the range of 0.4-0.7 fim. 

The host star is an M star with T^ff ~ 3800 K. There is a 
companion star 1.9" away with Teff ~ 3300 K. We infer a dis¬ 
tance to the system of 225 ± 50 pc, which, in turn, implies a 
projected physical separation of the two stars of ^430 AU. 
At this distance, the companion star may have been incidental 
to the formation of the planet, but perhaps was instrumental 
in driving it toward the host star via Kozai-Lidov cycles with 
tidal friction (Kozai 1962, Lidov 1962; Kiseleva et al. 1998; 
Fabrycky & Tremaine 2007). 

The minimum transit depth of < 0.14% sets an upper limit 
to the size of the underlying hard-body planet of 2.5 for 
an assumed radius of the host star of 0.57 ±0.067?©. This 
is consistent with the low surface gravity that is required to 
drive off metal vapors that could condense into dust (Rappa- 
port et al. 2012; Perez-Becker & Chiang 2013). In fact, it 
seems likely if the dust-emitting scenario we report for K2- 
22b is correct then Mars, Mercury, or even lunar sized bodies 
with surface gravities of 1/6 to 1/3 that of Earth are to be pre¬ 
ferred. 

We find that the dust tail of K2-22b has two properties that 
are distinct from those of KIC 125510 or KOI 2700b; (1) 
a leading dust tail (as opposed to a trailing one); and (2) a 
characteristic (e.g., exponential) scale length for the dust tail 
that is < 1/2/?*. The former requires dust transported to about 
twice the planet’s radius in the direction of the host star until 
it effectively overflows its Roche lobe, thereby going into an 
orbit that is faster than that of the planet. It is also necessary 
to have f3 (the ratio of radiation pressure forces to gravity) be 
< 0.02 which would be the case for very low luminosity host 
stars in combination with very small (< 0.1 ^m) or very large 
(> 1 ^m) dust particles (see Fig. fT4l) . The shorter dust tails 
could result from a combination of very low values of /3 and 
a short dust-grain sublimation time of < an hour. 

Mass loss rates in the form of high-Z material from this 
planet are likely to be 1.5 x 10^' g s“^ implying a lifetime 
of 20-70 Myr, for planet masses in the range of the moon’s 
to Mercury’s mass. 

The scenario described above for explaining the various 
features of the transit light curve and its variability is not en¬ 
tirely self-consistent. The requirement for particles to enter a 
leading tail is < 0.02 /rm which implies particle sizes of ei¬ 
ther a < 0.1 ^m or a > 1 pm. At the same time, the post-transit 
bump requires a forward scattering peak that is comparable to 
the angular size of the host star. This corresponds to particle 
sizes of ^ 1/2 pm. Finally, the color dependence of the tran¬ 
sits observed with the GTC implies a non-steep power-law 
size distribution with a maximum size of ^1/2 pm. Thus, 
larger particles could account for all three of these observa¬ 
tional pieces of evidence—^but only if a substantial fraction 
of the particles are large (e.g., ~ 1 pm). This requirement for 
such large particles seems somewhat unusual in comparison 

We note that Bochinski et al. (2015) were able to measure a small dif¬ 
ference in transit depths in KIC 1255b between g’ and z’ bands, and thereby 
concluded that the dusty effluents in this object contained a component of 
larger grains in the range of 0.25 - 1 /rm. 


TABLE 6 

Comparison of Dusty Planets 


Parameter (units) 

EPIC 201637* 

KIC 1255b2 

KOI 2700b3 

Arb (hr) 

9.146 

15.68 

21.84 

Depth (%) 

0-1.3 

0-1.4 

0.031-0.053 

Variability 

highly 

highly 

slowly 

d/R, 

3.3T0.2 

4.3 ±0.4 

5.9 ±0.4 

9* (deg) 

17 

13 

10 

Pre-bump 

weak 

yes 


Post-bump 

yes 

no 

no 

Stall (deg)'* 


~10-15 

~24 

Stall/S* 

;$o-5 

~0.77-l.l 

~2.4 

7i.o.t (K) 

3830 

4300 

4435 

Teq (K)5 

2100 

2100 

1850 

/3max 

0.05 

0.1 

0.07 

Mdust (10‘° g/s) 

20 

20 

0.15 


* Results from this work. 

^ Values from Rappaport et al. (2012); Brogi et al. (2012); Budaj 

(2013); van Werkhoven et al. (2014). 

^ Values from Rappaport et al. (2014). 

^ Approximate exponential tail length in degrees of orbital phase. 

with other known astrophysical collections of dust such as in 
the ISM (e.g., Mathis et al. 1977; Bierman & Flarwitt 1980), 
Solar-system comet tails (e.g., Kelley et al. 2013), lo (e.g., 
Kruger et al. 2003a; 2003b), and the Earth’s atmosphere (e.g., 
Liou 2002; Holben et al. 1998). 

The host star K2-22 shows photometric modulations with 
an amplitude of approximately 1% and a timescale of 15 days. 
It is interesting to note that the other two host stars to disin¬ 
tegrating planets also display photometric modulations at the 
several percent level, likely associated with star spots. It has 
been suggested that stellar activity may play an important role 
in modulating the process that generates the dust (Kawahara 
et al. 2013; Croll et al. 2015). The large spot modulations for 
the host stars of the candidate disintegrating planets is sugges¬ 
tive and warrants further investigation. 

Einally, we summarize in Table |6] some comparative prop¬ 
erties of the three known ‘disintegrating planets’. Hopefully, 
further patterns of similarity will emerge as more of these ob¬ 
jects are discovered. 

Discoveries of new disintegrating planets in the upcoming 
K2 fields could be potentially quite important. If found, these 
candidates are likely to orbit brighter stars than the host stars 
of the examples discovered to date. It would also be inter¬ 
esting to discover them orbiting a richer variety of host stars, 
especially since we have inferred from K2-22b that the low lu¬ 
minosity of the host star plays an important role in determin¬ 
ing the trajectories of the dust flowing from these very special 
planets. Euture follow-up observations with ground-based 
telescopes, both in photometric and spectroscopic modes, are 
likely to provide us with further information about the pro¬ 
cess that generates the dust emission, and give us better in¬ 
sights into the relevant physical processes in these extreme- 
environment systems. 
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